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Abstract 

We have identified a new type of transition from extended to localized vibrational states in one- 
dimensional graded elastic chains of coupled harmonic oscillators, in which the vibrating masses or nearest- 
coupling force constants vary linearly along the chain. We found that the derealization transition occurs 
at the maximum frequency of the corresponding homogeneous chain, which is in a continuous single band. 
Although each state in the localized phase, called gradon, can be regarded as an impurity localized mode, 
the localization profile is clearly distinct from usual impurity modes or the Anderson localized modes. We 
also argue how gradons may affect the macroscopic properties of graded systems. Our results can provide 
insights into many analogous systems with graded characters. 

PACS numbers: 63.20.Pw, 63.20.Dj, 63.22.+m 
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I. INTRODUCTION 



Vibrational normal modes and their localization properties in inhomogeneous or disordered 
media have attracted much interest . 'i 2 i 3 i 4 i 5 i 6 Notable examples include disorder-induced local- 
ized states (i.e., the Anderson localization of lattice vibrations), fractal-induced localization (i.e., 
fractons 7 ), modes contributing the boson peak, 8 the connection of excess of low-frequency modes 
with the jamming-unjamming transition in particle packings, 9 intrinsic localized modes (ILMs) 
caused by nonlinearity and/or discreteness in systems, 10 and the well-known defect modes in 
crystals including impurities^ etc. Surface phonon modes in semiconductors are also typical 
localized vibrational excitations in which atomic displacements decrease exponentially into the 
bulk. 12 These localized vibrational modes strongly affect properties of condensed matter. It is cru- 
cial to reveal characteristics of localized modes for understanding thermal, elastic, and acoustic 
properties and/or electron-phonon interactions in materials. 

From the viewpoint of localization mechanisms, there exist two types of localized modes in 
noninteracting systems. The first type is a consequence of interference of coherent vibrational 
waves due to diffusive scattering. Localized modes of this type include the Anderson localiza- 
tion of classical wave o 13 ' 14 and fractons in fractal systems. 7 In these cases, localization occurs 
whenever the elastic mean free path becomes comparable to the inverse wave vector of the wave 
(the Ioffe-Regel condition). In particular, the Anderson localization of elastic waves in disordered 
solids has been extensively studied theoretically and experimentally. 1 ^ 16 ^ 7 Now, we have obtained 
vast amount of knowledge on the Anderson localization and the Anderson transition, which in- 
volve similarities to and differences from quantum particle localization. 18 The mode pattern of an 
Anderson localized state possesses an exponential tail and the position of the localization center 
is unpredictable. The scaling theory 19 shows all vibrational modes are localized in one- and two- 
dimensional disordered systems while one has the localization-delocalization transition in three 
dimensions, which is the same with the Anderson transition of quantum particles in the orthogonal 
universality class. 20 The derealization transition occurs in a continuous band without any singu- 
larity of the spectrum at the transition frequency uj c . The localization length £ diverges as a power 
law of \u — u c \. The other type of localizations is due to confinement by impurities such as defect 
modes, surface modes, or ILMs. 21 This type of localized modes is rather trapped by local poten- 
tials. The mode pattern of such an impurity localized mode has also an exponential tail and the 
position of the localization center is located at the impurity site. For instance, the light-mass impu- 
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rity mode of a one-dimensional elastic chain has the mode pattern of u n oc [—m / (2m — m )]' n ', 
where m and m (< m) are the host and impurity masses, respectively, and n is the site number 
(n = represents the impurity site). The frequency of this light-mass impurity mode is given 
by u>? = m 2 dj^ ax /m (2m — m ), where u max is the maximum frequency of the homogeneous 
chain with mass m for all sites. Since m < m, u- imp is always larger than u max , which im- 
plies that the light-mass impurity localized mode has an isolated spectrum outside of the band. 
This is considered as a common feature of impurity localized modes. Therefore, it seems that no 
localization-delocalization transition induced by impurities occurs in a single band. This is, how- 
ever, not obvious. If a derealization transition in a continuous spectrum is realized by impurity 
modes, properties of localized vibrations and the transition would be quite different from those for 
conventional localized modes and transitions. 

In this work, we consider vibrational problems of a graded elastic systems which are modelled 
by graded spring-mass networks. Studying one-dimensional (ID) systems with graded masses 
and graded force constants, we show a new type of transition from extended to impurity localized 
modes. Some characteristic properties of these localized modes are also elucidated. Graded elastic 
systems can be found everywhere in nature. Mass density of an elastic body subject to a gravi- 
tational field is low near the top of the object but high near the bottom, such as sediment at the 
bottom of ocean, which can be regarded as an example of graded elastic systems. Similar localized 
modes to our results can be expected in it. A carbon nanotube which can be regarded as folding 
a graphite sheet into a tube is another example of graded elastic materials. 22 The force constant 
can be stronger at the central part of the tube but weaker at both ends. Thus there is a gradient of 
force constant along the nanotube. Particularly, by virtue of the analogy of the vibrational problem 
with many other problems, such as electromagnetic wave propagation 2 ^! in which the index of 
refraction plays the same role as the mass, one would expect similar localizations in other type of 
graded structures such as quantum or photonic systems. It is possible to apply our model to graded 
plasmonic nanoparticle chains which may find prominent applications in nanooptics. In fact, there 
have been many efforts directed to graded systems in both artificial composites 25 and materials 
in nature (e.g., biological cells 26 ). One further notices that in contrast to a simple sharp inter- 
face, it has been studied that the presence of gradient in a interface leads to frequency-dependent 
properties; Although we treated ID chains in this work, our findings form the basis for higher 
dimensional networks in which results are more complicated and intriguing and we will present 
them elsewhere. 29 
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In what follows, we first briefly describe, in Sec. II, the graded elastic network model, and 
establish the formalism for studying eigenmodes. In Sec. [nil we present our results mainly in the 
graded mass model and similar results are in the graded force constant model . Finally conclusions 
and discussions on our results are given. 



II. MODEL AND FORMALISM 



The elastic vibrational problem of a uniform three-dimensional system is simply modelled by 
a cubic lattice of harmonic oscillators where different components of the displacement decouple. 8 
However, when there exists a uniaxial gradient of coupling constants (or masses), the equation 
of motion is only partially separable. Assuming that the network is periodic in the directions 
perpendicular to the gradient direction (e.g., z-axis), one reaches an effective one-dimensional 
problem. 29 Without loss of generality, let us start with a one-dimensional spring-mass model as 
depicted in Fig. [lj The mass M n and the spring constant K n connecting to the right nearest- 
neighbor are labelled sequently with n = 1, 2, • ■ ■ , N. The motion of each mass is governed by 
the following equations 1 

M n u n = K n _i(u n _i - u n ) + K n (u n+1 - u n ), (1) 

where u n denotes the scalar displacement of mass M n at the site n, and the derivative u n is taken 
with respect to time t. We omit damping and further couplings to avoid complications, the in- 
clusion of which do not change the results qualitatively. Decomposing the amplitude u n by the 
eigenvectors e(A) of the mode A as u n = ^ x Q\{t)e n (\), where Q\(t) is a time-dependent ex- 
pansion coefficient which behaves as Q\(t) ~ exp(— iuj\t), we simplify Eq. (HJ to the eigensystem 
of 

-uj 2 e = T e, (2) 

where the matrix T characterizes the mass-weighted coupling interactions of the neighboring 
masses connected by the springs that form a lattice with lattice constant a = 1. The system 
size L is then given by aN. Assuming free boundary conditions, the N x N matrix T becomes 
tridiagonal and with elements 



T 



'1 — S n l)K n -l , (1 — 5 n N)K n 



K n 

&mn j^j- (^mn-1 ^mn+l)i 



(3) 



where 5 is the Kronecker delta. In this regard, further incorporating the next-nearest coupling force 
constant K" ext results a quintuple matrix T. We have also attempted periodic boundary conditions, 
i.e., let M n connected to M t by K n , very similar results are observed. The nontrivial N solutions 
of e and lu 2 in Eq. © are thus determined by diagonalizing T. 

In this paper, we will consider two different cases in a one-dimensional network: 

(1) The graded mass case with 

M n = M - C M n/N, (4) 

and fixed K [see Fig.[TJa)]. In this case, the node mass decreases linearly from the left end site 
n = 1 to n = N. 

(2) The graded force constant case with 

K n = K + C K n/N (5) 

and fixed M [see Fig. E^b)]. Here, M and K are the maximum mass and the minimum force 
constant, and C M and C K are the gradient coefficients in the two models, respectively. Note that 
the conditions of Cm = and Ck = recover the homogeneous chain. It is also emphasized 
that our graded models can be easily extended to infinite N (namely, infinite system size L). The 
infinite graded chain does not correspond to the continuum graded chain, but to the infinitesimal 
gradient limit. 

HI. NUMERICAL RESULTS 

In this section, we mainly concentrate on the graded mass model with iV = 1000. A further 
increase in N does not change our results (except for Fig. |5] as shown later). The dispersion 
relation of the vibrational modes in a crystal gives rich information about properties of phonons. 
In the graded elastic network, however, the wavenumber is not conserved because of the lack of 
translational symmetry in the system. This does not allow us to consider the dispersion relation 
for our model. Nevertheless, the relation between the eigenfrequency uj\ and the mode index A in 
ascending order with respect to u\, which can be regarded as a pseudo-dispersion relation, helps 
us to interpret vibrational modes in the graded system. Figure El shows the eigenfrequency lo\ 
as a function of the eigenmode index A for gradient coefficients Cm = 0.0, 0.2, 0.5, and 0.7. 
M = 1, Kq = 1 and Ck = are used. Note that the relation for the homogeneous elastic chain 
(C M = 0.0) represents the usual dispersion relation. It is clear that the presence of the gradient 
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alters the (pseudo-) dispersion relation, particularly beyond the frequency uj c = 2^K$/Mq = 2.0 
(dashed horizontal line) which is the maximum eigenfrequency of the homogeneous elastic chain 
with M n = M for any n (the solid curve). The number iV g of high-frequency excitations with 
uj > uj c depends on the gradient coefficient C M as shown by the insert of Fig. El The maximum 
frequency of the graded chain is given by u max = 2^K /M N (2^jK N /M ) in the graded mass 
(force constant) model. Let us call these high energy excitations "gradons" which are clearly 
distinct from low frequency excitations as we will show later. The fact that uj c is a characteristic 
frequency in the graded elastic chain is more explicitly found in the density of states (DOS) defined 
by 

A 

The DOS shown in Fig. |3] clearly exhibits a singular maximum peak at uj c = 2.0 for C M > 0. This 
behavior persists for various gradient coefficients Cm- The low-frequency DOS is in accordance 
with the usual asymptotic Debye density of vibrational states, both in the homogeneous case and 
the graded case. From the DOS, we expect that modes with uj > uj c are qualitatively different from 
modes with u < u c . In order to clarify this point, we plot the mode patterns e(A) in Fig.HJ where 
three typical modes with u = 2.38223(> uj c ), 1.99845(« uj c ), and 0.39309(< uj c ) are presented. 
These mode patterns show that modes with uj > uj c (namely, gradons) are confined at the "light" 
mass side, whereas those with uj < uj c (phonons) are extended. As will be shown later, these con- 
fined modes are qualitatively different from conventional localized modes. Nevertheless, we use 
the word localized modes for these confined states, because the present model exhibits a drastic 
change in spatial extent of modes. One can find a beating pattern in the wave profile [see Fig.Etb)], 
which implies that the mode is composed of plural wavenumber components. The degree of mix- 
ing of different wavenumber components can be quantified by the full width at half maximum 
(FWHM) of the Fourier amplitude \Fj(e)\ = |tt J2n=i e « ex P \ 2j[i U ~ - 1)/N]/VN\. In- 
deed, the Fourier transform of the mode pattern in Fig.Ua) has a finite width around the average 
wavenumber as shown in Fig. 5(a). Figure 5(b) gives the FWHM as a function of uj for several val- 
ues of Cm- The FWHM shows singularity at uj = uj c as in the case of the DOS. For well-extended 
states for uj uj c , the FWHM must be very small, which implies that we can approximately define 
the Bloch wavenumber k for modes in this very low frequency regime. It is then possible to regard 
Fig.|2]as a usual dispersion relation between uj and k for uj uj c . 

In order to quantify the degree of localization of excited modes, we consider the inverse partic- 
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30.31 



(IPR), 



N 



P 



-i 



n=l 



(7) 




The IPR defined by Eq. © gives a measure of the number of sites contributing significantly to 
a given eigenstate. If the state is strongly localized, becomes proportional to l/£ d , where 
£ is the localization length, while for a fairly extended state it is proportional to 1/L d because 
all sites contribute to the state. We show in Fig. 6 the IPR of vibrational modes in the graded 
mass model as a function of the eigenfrequency for various gradient coefficients Cm- At low 
frequencies {uj < uj c ), small and almost flat IPR indicates fairly well-extended excitations across 
the whole network. The IPR increases dramatically, however, as the frequency exceeds uj c = 2.0. 
This shows that modes with uj > uj c are spatially localized. These are in accord with the excitation 
shapes depicted in Fig. |4j As seen in the inset of Fig. 6, the small step of the IPR at uj = uj c is 
somewhat broad. This is a consequence of the finite size effect of our N = 1000 system. We 
confirmed that the step becomes sharper as increasing N (figure not shown here). 

Due to the duality between force constant and mass, the graded force constant model repre- 
sented by Eq. © is essentially the same with the graded mass model we discussed above. Indeed, 
we basically observed the same results as those in the graded mass model. In Fig|7J for exam- 
ple, we present three typical excitation states with u = 2.27844(> u c ), 1.99911(~ uj c ), and 
0.3776 18 (< uj c ) for the graded force constant model with C K = 1.0. These mode patterns are 
quite similar to those in the graded mass model shown in Fig. |4] Further results on the pseudo- 
dispersion relation, the DOS, the FWHM, and the IPR show also similar behaviors to those in the 
graded mass case. 

From all the results presented above, we can conclude that an infinite graded elastic chain 
exhibits the localization-delocalization transition at the frequency uj c being the maximum eigen- 
frequency of the homogeneous elastic chain with M n = M (or K n = K ) for any n. In the 
frequency regime uj < uj c , we have phonon-type modes extended over the whole system, while 
modes with to > uj c (namely, gradons) are localized in the lighter mass side (or harder force con- 
stant side). It should be noted that this derealization transition (the phonon-gradon transition) 
occurs in a single band as shown in Fig. 3. The mechanism of gradon localization is considered 
as follows. A mode with a frequency uj larger than cu c cannot have amplitudes in the heavier (or 
softer) region because masses (or force constants) in this region are too heavy (or too soft) to 
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vibrate with the frequency w. As a consequence, the mode pattern has finite amplitudes only in 
the lighter (or harder) region, while no amplitude in the heavier (or softer) region. The boundary 



The fact that values of n c of the modes shown by Fig. Ufa) and Fig. 7(a) are predicted by Eq. © 
as n c = 423 for Fig.Uta) and 298 for Fig. 7(a) shows the validity of our interpretation of gradon 
localization. This mechanism of localization is essentially the same with that of impurity localized 
modes. Therefore, gradons belong to a kind of confined modes by impurities, and the phonon- 
gradon transition in a single band is induced by impurity localized modes. 

Properties of this new type of the derealization transition are largely different from usual tran- 
sitions occurring in a single band, such as the Anderson transition. First of all, the gradon transition 
can occur even in a one-dimensional graded system, while the Anderson transition requires three 
or more dimensions. This is not surprising, because the scaling theory cannot be applied to the 
gradon transition. Furthermore, a localized mode in a usual sense has a finite localization length 
£ even if the system size is infinite. On the contrary, the localization length of a gradon, which is 
defined as the size of vibrating region and approximately given by N — n c , becomes infinite if the 
system size N tends to infinity as seen from Eq. ©. Gradon modes are localized only in the sense 
that a part of the whole system has vibrational amplitudes. The ratio of the vibrating region to the 
whole system is given by 1 — M (l — 7 2 )/C M for the graded mass model and 1 — Kq{^~ 2 — 1)/Ck 
for the graded force constant model, where 7 = uj/uj c . Although it is known that there exist one- 
dimensional derealization transitions in correlated disordered systems, 32 the gradon transition is 
distinguished from this type of transitions at this point. The fact that the DOS shows the singu- 
larity at the transition frequency is also a feature of the phonon-gradon transition, which contrasts 
strikingly with the Anderson transition without any spectral singularity at the critical point. For 
the graded mass model, each state in the vibrational DOS above co c corresponds to a light-mass 
impurity mode mentioned in Sec. I. The linear distribution of light-mass impurities in the graded 
mass model forms a continuous spectrum of impurity modes above co c . The whole spectral struc- 
ture shown in Fig. 3 is constructed by connecting two different spectra, the phonon-type excitation 



between these two regions is determined by uo = 2^J K n /M n . Therefore, the localization center 
site n c which is the position with the maximum amplitude [see Fig. Eta)] is given by 



n. 




graded mass model, 



graded force constant model. 



(8) 



8 



spectrum and the impurity mode one above uj c . 

As seen from Figs. |4] and |7J mode pattern profiles of gradon excitations are highly asymmetric 
due to the asymmetricity of the graded elastic network. A gradon mode in the heavier (or softer) 
side of the localization center n c has an exponential tail with a large damping factor d^ 1 as ~ 
exp[— (n c — n)a/d ). The width of the localization front (i.e., the decay length) d seems to be 
independent of uj. On the contrary, in the lighter (or harder) side of n c , vibrational amplitudes 
decrease very slowly (almost constant) as increasing n — n c . We should remark that the gradon 
transition remains sharp even for a finite system. This is because the width d is much smaller than 
the system size L. This also contrasts to a conventional transition described by the scaling theory, 
in which the finite system size reduces the transition to a broad crossover. 

IV. CONCLUSION AND DISCUSSION 

In summary, we identified a type of transition from extended to localized states in the one- 
dimensional graded elastic network described by Eq. (HJ with Eq. © or ©. The derealization 
transition occurs at uj = uj c , where uj c is the maximum frequency of the homogeneous network 
with M n = M (or K n = K ) for any n. The transition is totally distinguished from conven- 
tional disorder induced derealization transition such as the Anderson transition in the following 
sense: (1) The derealization transition occurs even in a one-dimensional system in contrast to the 
Anderson transition requiring three or more dimensions. (2) Localized vibrational states which is 
called gradons in the graded elastic network are classified into impurity localized modes though 
the spectral density of states (DOS) of gradon modes is continuous. (3) The localization length 
of a gradon becomes infinite if the system size tends to infinity. (4) At the transition frequency 
uj c , the DOS shows a singularity, which is a consequence of the connection of the phonon-type 
excitation spectrum to the one of impurity mode at uj = uj c . The DOS of infinite graded chain can 
in fact be expressed as an average of the DOS of homogeneous chains over a range of masses. (5) 
The mode pattern of a gradon excitation is highly asymmetric around the localization center which 
is the position with the maximum amplitude. (6) The gradon transition remains sharp even for a 
finite system because the width of the localization front is much smaller than the system size. 

It would be of great interests to study the general graded profiles. In this regard, further calcu- 
lations with other graded profiles have been attempted and qualitatively similar results have been 
obtained. For instance, the power-law graded profile still confirms the existence of the transition 
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from phonons to gradons. However, the transition frequency uo c ^ 2 any longer; it depends on 
the detailed graded profile. Higher-dimensional graded networks are more realistic models. These 
networks are possibly in two or three dimensional crystals, in which the vibrational modes can be 
partitioned in the dispersion relations, due to the presence of bounded surface elastic waves.— In 
current considerations, we have not included any damping or non-neighboring couplings and open 
boundary conditions are enforced. We believe the inclusion of further coupling or damping will 
not change the conclusions much, although inclusion of a damping can change the normal modes 
with an overall time decaying factor. In fact, further calculations with the next-nearest coupling 
K" ext in the model show enhancement of the degree of localization, increasing of the transition 
frequency uj c but decreasing of the number of localized modes. 29 The localization tail still ex- 
hibits an exponential rather than power-law decay, the latter of which occurs only in systems with 
long-range couplings. 29 Localized states in high dimensions with vector displacements are also 
possible, but the bending and stretching force constants can introduce a new characteristic length 
scale and more complicated spectrum in such systems. 

It is also natural to think about how physical quantities are affected by the existence of gradons. 
Since the DOS possesses a sharp peak at the phonon-gradon transition frequency uj c , we could 
expect a remarkable structure near T = T = fiw c / K B in the heat capacity, just as the hump of the 
specific heat capacity due to boson peak. 8 Furthermore, we believe the temperature dependence 
of the thermal conductivity k(T) is strongly affected by gradons due to the fact that the diffusion 
constant for localized gradon modes is zero, which facilitates locking of the heat flux. It is because 
if we neglect the anharmonic interaction between phonons and gradons, the thermal conductivity 
would become constant for T ^> T . One cannot, however, rule out the anharmonic interactions 
between phonons and gradons at high temperatures, and the decay processes from gradons to 
phonons give an additional thermal conductivity. As a result, the temperature dependence of k{T) 
for T ^> T would be different from n(T) for crystals, which provides experimental evidence 
of gradon modes in the graded elastic networks. Further incorporating the gradon results with 
diatomic model, anharmonicity of nearest neighbor interaction, and rotational degree of freedom 29 
is also rewarding. We expect experiments be done to observe gradons, or even its analogous like 
spin-wave gradons, quantum gradons, and photonic gradons etc. The present results and those 
in higher dimensions 29 will have implications in surface acoustic wave applications as diverse as 
touchscreens, pizeoelectric material design, and integrated microresonator systems etc, and also 
provide benefits to a wide spectrum of problems such as earthquake study or oil exploration. 
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Figure Captions 

FIG. 1. Schematic pictures of the graded elastic networks, (a) The graded mass model, and (b) the 
graded force constant model. 

FIG. 2. The pseudo-dispersion relation for the graded mass models with gradient coefficients 
C M = 0.0, 0.2, 0.5, and 0.7 and M = 1.0. C M = corresponds to the homogeneous case [i.e., 
k = nix / L and uo = 2uo sin(fca/2)]. Inset shows the number of gradons iVg as a function of C M - 

FIG. 3. The densities of states for the graded mass model with M = 1.0 and various C M . 

FIG. 4. Three typical excitation states of the graded mass model with M = 1.0 and C M = 0.7 at 

(a) uo = 2.38223(> u e ), (b) uo = 1.99845(« uo c ), and (c) uo = 0.39309(< uo c ). 

FIG. 5. (a) The Fourier transform of the mode pattern shown by Fig. 4(a). (b) The full width at 
half maximum of Fourier transforms of mode patterns of the graded mass model with M = 1.0 
as a function of uo. 

FIG. 6. The inverse participation ratio (IPR) of modes excited in the graded mass model with 
M = 1.0 plotted as functions of the eigenfrequency. Inset shows the magnified IPR around 
uo = 2.0. 

FIG. 7. Three typical excitation states of the graded force constant model with K = 1.0 and 

C K = 1.0 at (a) uo = 2.27844(> uo c ), (b) uo = 1.99911(« uo c ), and (c) uo = 0.377618(< uo c ). 
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Eigenmode index 



FIG. 2: /Xiao, Yakubo, and Yu 
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FIG. 3: /Xiao, Yakubo, and Yu 
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